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Abstract 

Using 14 million ip(2S) events collected with the Beijing Spectrometer (BESII) at the Beijing 
Electron-Positron Collider, the branching fractions of ip(2S) decays to pfnr~ and pmr + and the 
branching fractions of the main background channels ip(2S) — > pnir~ir® , ip{2S) — > 7XcO — * ~ipnir~ , 
tp(2S) — ► 7Xc2 — ► 7pn7r _ , and ip(2S) — > 7Xcj — ► 7pn7r _ 7r° are determined. The contributions of 
the iV* resonances in ip(2S) — > pn7r _ + c.c. are also discussed. 

PACS numbers: 14.20. Gk, 13.75.Gx, 13.25.Gv 
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I. INTRODUCTION 



From perturbative QCD (pQCD), it is expected that both J/0 and -0(25) decaying 
into light hadrons are dominated by the annihilation of cc into three gluons, with widths 
proportional to the square of the wave function at the origin |\I/(0)| 2 [1]. This yields the 
pQCD "12% rule" 

Q h = « ^ (25) - e+e ' « 12%. 

Bj/^-th Bj/^ e + e - 

The violation of this rule was first observed in the pn and K* + K~ + cc. decay modes 
by Mark-II [2]. Following the scenario proposed in Ref. [3], that the small i>{2S) — > pit 
branching fraction is due to the cancellation of the 5- and D-w&ve matrix elements in 0>(25) 
decays, it was suggested that all 0(25) decay channels should be affected by the same 5- and 
.D-wave mixing scheme, and thus all ratios of branching fractions of 0(25) and J/tp decays 
into the same final state could have values different from 12%, expected between pure 15 
and 25' states [4]. The mixing scenario also predicts -0(3770) decay branching fractions since 
the -0(3770) is a mixture of S- and D-wave charmonia, as well. Many channels of J/-0, 
■0(25), and 0(3770) decays should be measured to test this scenario. 

A very important source of information on nucleon internal structure is the N* mass spec- 
trum, including production and decay rates. Because of its importance for the understanding 
of nonperturbative QCD, a series of experiments on N* physics with electromagnetic probes 
(real photons and electrons with space-like virtual photons) are being performed at facilities 
such as JLAB, ELS A at Bonn, GRAAL at Grenoble, and SPRING8 at JASRI [5]. They 
have already produced some results [6, 7]. However, our knowledge on N* resonances is still 
poor. Even for the well-established lowest excited state, the A/"* (1440), properties such as 
mass, width, and decay branching fractions still have large experimental uncertainties [8]. 
Another outstanding problem is that, in many of its forms, the quark model predicts a 
substantial number of N* states around 2 GeV/c 2 , which have not yet been observed [9]. 

Recent studies of N* resonances have been performed using J ftp events collected at the 
Beijing Electron-Positron Collider (BEPC) [10, 11], providing a new method for probing 
this physics, and a new N* peak with a mass at around 2065 MeV/c 2 was observed [11]. 
This may be one of the "missing" N* states around 2 GeV/c 2 . However, due to its large 
mass, the production of this A^*(2065) in J/0 decays is rather limited in phase space. A 
similar search for it in 0(25') — > ppn has been performed [12], where there is a faint but not 
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statistically significant accumulation of events in the pn invariant mass spectrum at around 
2065 MeV/c 2 . 

In this paper, we study ip{2S) — > pun' + c.c. and their main background channels, 
determine branching fractions, test the 12% rule, and study N* resonances in the Ntt system. 

II. BES DETECTOR AND THE DATA SAMPLE 

BESII is a large solid-angle magnetic spectrometer which is described in detail in Ref. [13]. 
The momentum of charged particles is determined by a 40-layer cylindrical main drift cham- 
ber (MDC) which has a momentum resolution of <t p /p=1.7%a/1 + p 2 (p in GeV/c). Particle 
identification (PID) is accomplished using specific ionization (dE/dx) measurements in the 
drift chamber and time-of-flight (TOF) information in a barrel-like array of 48 scintillation 
counters. The dE/dx resolution is <JdE/dx — 8.0%; the TOF resolution for Bhabha events 
is cttof = 180 ps. Radially outside of the time-of-flight counters is a 12-radiation- length 
barrel shower counter (BSC) comprised of gas tubes interleaved with lead sheets. The BSC 
measures the energy and direction of photons with resolutions of cr E /E ~ 21%/\/E (E in 
GeV), <j0 = 7.9 mrad, and a z = 2.3 cm. The iron flux return of the magnet is instrumented 
with three double layers of proportional counters (MUC) that are used to identify muons. 

In the analysis, a GEANT3-based Monte Carlo (MC) simulation program (SIMBES) with 
detailed consideration of the detector performance is used. The consistency between data 
and Monte Carlo has been checked in many high purity physics channels, and the agreement 
is reasonable [14]. For the MC generators, the angular distribution for tp(2S) — > 7% c j is 
simulated assuming a pure El transition, and uniform phase space is used for the other 
decays. 

The data sample used for this analysis consists of (14.0 ± 0.6) x 10 6 ip(2S) events taken 
at y/s = 3.686 GeV [15]. Backgrounds are estimated using an inclusive ip{2S) decay MC 
sample generated by LUNDCRM [16] with the same size as the ip(2S) data. 

III. EVENT SELECTION 

For the signal channel i/j(2S) — > pn-n~ +c.c. and the background channels, we reconstruct 
two charged tracks, and the neutron and antineutron are not measured. However, since 
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most antineutrons annihilate in the detector (mainly in the BSC) and most neutrons pass 
through the detector without interaction, these signatures are used to suppress backgrounds 
by requiring a neutral cluster in the expected antineutron direction and requiring nothing 
in the neutron direction. 

A neutral cluster is considered to be a good photon candidate if the following requirements 
are satisfied: it is located within the BSC fiducial region; the energy deposited in the BSC is 
greater than 50 MeV; the first hit appears in the first 6 radiation lengths; the angle between 
the cluster development direction in the BSC and the photon emission direction from the 
beam interaction point (IP) is less than 37°; and the angle between the cluster and the 
nearest charged particle is greater than 15°. 

Each charged track is required to be well fit by a three dimensional helix, to originate 
from the interaction point, V xy = ^JV£ + Vy < 1 cm, \V Z \ < 15 cm, and to have a polar 
angle | cos#| < 0.8. Here V x , V y , and V z are the x, y, and z coordinates of the point of closest 
approach to the beam axis. The TOF and dE/dx measurements for each charged track are 
used to calculate Xpid(^) values and the corresponding confidence levels Probpwii) for the 
hypotheses that a track is a pion, kaon, or proton, where i (i — ir/K/p) is the particle type. 

We use the following common selection criteria for all channels: 

1. The number of charged tracks in the MDC is two with net charge zero, and the 
difference in V z between the positive and negative charged tracks |Vjt — V~\ is required 
to be less than 3 cm. 

2. For each charged track, the particle identification confidence level for a candidate 
particle assignment is required to be greater than 0.01. For a proton we also require 
Probpioip) > Probpioin) and Probpinip) > P r obpio(K), and for a n, we require 
Probpi D (n) > Probpinip) an d Probpio^) > Prob PID (K). 

3. The energy of the positive charged track observed in the BSC should be less than 0.7 
GeV in order to eliminate ip{2S) — > e + e~ and Bhabha events. 

4. M p7r -(or p7T+) > 115 GeV/c 2 to remove background channels with A or A. 

For ip(2S) — > pnn~ + c.c. when there are photon candidates in the events, a < 10° is 
required, where a is the angle between the nearest neutral cluster and the missing momentum 
direction of the two charged tracks. When there are no photon candidates in an event, the 
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value of a is set to zero degree. After all above requirements, we obtain the pm~ and p7i + 
missing mass distributions shown in Fig. 1, where there are clear signals of n and n at around 
0.94 GeV/c 2 . The background level increases from about 10 events per 20 MeV/c 2 bin at 
0.6 GeV/c 2 to about 30 events at about 1.2 GeV/c 2 . 
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FIG. 1: Missing mass distributions of (a) pir~ and (b) pir + of ip(2S) — > pfnr~ +c.c. candidate events. 
The squares with error bars are data, the blank histograms are the fit, the shaded histograms 
are normalized backgrounds from ip(2S) — ► 7XcJ,XcJ — > pnTT~ + c.c, tp(2S) — > pn7t~7T° + c.c., 
ip(2S) — > 7Xcj, Xcj — ► pn7r _ 7r° + c.c, and the dashed curves are background shapes from the fit. 



IV. BACKGROUND ANALYSIS 

Backgrounds are studied using the ip{2S) inclusive decay MC sample. By applying the 
same selection criteria, it is found that for ip(2S) — > pfnr~ the main backgrounds are from 
ip(2S) -> iXc.j,Xcj -> pnic'ii , 4){2S) -> pnn-71 , and ^(25) -> 7XcJ, XcJ -> P^vr". All 
these modes have not been measured previously. Here we measure these channels. 

For these three background channels, in addition to the above selection criteria, we also 
require: 

1. For ip(2S) — > ^XcjiXcj pn7c~7[° events, we do a one-constraint (1C) kinematic fit 
to '0(25') — ► 777pn7r~ looping over all photon candidates. The combination with the 
minimum x 2 is selected, and Prob\c > 0.01 is required. After the 1C fit, the direction 
of the n is determined. Similar requirements are imposed on ip(2S) — > / ~f / ~fpn7c~ and 
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i/j(2S) — > 7pn7r final states in selecting ip(2S) — > pn7r tt° and i/j(2S) — > 7XcJ, Xcj — * 

2. The angle between the direction of n and one of the neutral clusters not already used 
in the kinematic fit should be less than ten degrees. 

A. ip(2S) —> JXcJ, Xcj — > pniv-n 

For this channel, we require four photon candidates (three of them from the radiative and 
7T° decays, the other from the interaction between the antineutron and the detector material) 
and the angle between the two photons from tt° decays be greater than eight degrees to 
remove the background from split-off fake photons. In order to suppress the background 
from ip(2S) — > 7r°7r°J/ip : J/t/j — > pnn~, the requirement \m pfi7r ~ — 3.097| > 0.05 GeV/c 2 is 
used. After applying the above selection criteria, the 77 invariant mass spectrum, shown in 
Fig. 2, is obtained, and |m 77 — 0.135| < 0.03 GeV/c 2 is required to select ir° candidates. 



FIG. 2: The 77 invariant mass distribution for selected ip(2S) —> jpnit 77 candidate events. 



date events is shown in Fig. 3. The shaded histogram shows the main backgrounds from 
ip{2S) — > 7i 7T°J/ilj,J/ip — > pnir~ and ip(2S) — > pfnr~n°, which have been normalized to 
data using the branching fractions of J/ip — > pnir~ and ip{2S) — > pnn'TT from Ref. [11] 
and our measurements (see Sec. IV B). From Fig. 3, no clear Xcj signals are seen. Since the 
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background in this channel is very complicated, we set an upper limit for Xcj pnir~ir by 
subtracting the known backgrounds from the total number of observed events and obtain 
45 ± 9 events for pfnr^ir invariant mass greater than 3.2 GeV/c 2 . Assuming a Gaussian 
distribution, the upper limit at the 90% C. L. is 57. Using the MC simulated efficiency e x 
of (4.00 ± 0.09)%, we obtain: 

2 jyup 

Y^BmS) - 1Xc j,Xcj -pmr-vr ) < = 1.2 xlO" 4 , 

~ N^2S) ■ ei • h ■ h ■ (1 - s) ■ -> 77) 

at the 90% C. L. Here f\ and are efficiency correction factors (see Sec. V, items 6 and 8) 
and s is the systematic error (see Table I). 
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FIG. 3: The pnTT~^^ invariant mass distribution for selected ip(2S) — > 'ypnn^-n candidate events. 
The shaded histogram is the sum of ip(2S) — ► 7r°Tr° J/ ip, J/ tp — > pnir~ and tp(2S) — > pnn^-n 
backgrounds which have been normalized to data according to their branching fractions. 



B. ip(2S) — > pmr-7r 

For this channel, we require the number of selected photon candidates equals three (two 
of them come from 7r° decays and the other from the interaction of the antineutron with the 
detector material) and the angle between the two photons from 7r° decays be greater than 
eight degrees to remove the background from split-off fake photons. The 77 invariant mass 
distribution after applying all the above selection criteria is shown in Fig. 4, where a clear tt° 
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signal can be seen. From the exclusive MC simulation, it is seen that possible backgrounds, 
except ip(2S) — > 7Xcj, XcJ pn^'n , have no peak at the n° mass in the 77 invariant mass 
distributions, and therefore they will not contribute to the number of n° events in fitting 
the 77 invariant mass distribution. For ip(2S) — > ^Xcj,Xcj pfnr~7r° background, after 
normalizing to data according to the branching fraction measured above, the contribution 
of this background is so small that it can be neglected. 

By fitting the 77 invariant mass spectrum with a MC determined signal shape and a 
2nd order Legendre polynomial for background, as shown in Fig. 4, 135 ± 21 ip(2S) — > 
pfm~n° candidate events are obtained. The statistical significance is 8.1a, and the detection 
efficiency e 2 for this decay mode is (3.21 ± 0.08)%. We obtain: 



N sig 



(3.18 ±0.50) x 10" 



N^s) ■ €2 • fi • fs • B (n° -> 77) 
where fi, f% are the efficiency correction factors, and the error is statistical. 




0.10 0.20 0.30 0.40 
yy invariant mass (GeV/c 2 ) 



FIG. 4: The 77 invariant mass distribution of the selected ip(2S) — > jjpnir candidate events. 
The squares with error bars are data, and the histogram is the fit. 



C. tp(2S) — > 7Xcj, XcJ — >• pnn 

For this channel, we require two photon candidates (one comes from radiative decay 
and the other one from the interaction between the antineutron and the detector material). 
After applying all the above selection criteria, the pnn~ invariant mass distribution, shown in 
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Fig. 5, is obtained. Here Xco and Xc2 signals can be seen, but the Xci signal is less significant. 
In fitting, we neglect the Xci signal, and 85 ± 18 and 80 ± 16 events are obtained by fitting 
the pnn~ invariant mass spectrum with MC simulated Xco and Xc2 signal histograms and a 
background shape. Here the background shape includes a 2nd order Legendre function and 
a normalized histogram obtained from background channel ip(2S) — > pnrc~ measured in this 
work. 

The statistical significance of Xco and Xc2 are both 4.2a, and the detection efficiencies for 
Xco (^3) and Xc2 (q) are (5.77 ±0.11)% and (6.14 ± 0.11)%, respectively. We obtain: 

N sig 

B^(2S) - 7Xeo, Xco - pnn~) = - x " n \ - = (1.10 ± 0.24) x 10~ 4 , 

^4>(2S) • £3 • Jl • J3 
N sig 

B(ij(2S) - 7Xe2, Xc2 - pnn-) = x ^ p ™ = (0.97 ± 0.20) x 10" 4 , 

^4>(2S) • £4 • Ji • h 

where fi and / 3 are the efficiency correction factors, and the errors are statistical. 

30 1 1 




FIG. 5: The pnir~ invariant mass distribution of the selected ip(2S) — > ^pnn~ candidate events. 
The squares with error bars are data, the histogram is the fit, and the dashed curve is the back- 
ground shape from the fit. 

Since the neutron passes through the detector without interaction in almost all cases, the 
background level will be higher and more complicated than for the antineutron final state, 
so we do not try to measure the charge conjugate modes of the above three channels. In the 
following analysis, we assume the branching fractions for these three channels equal to their 
charge conjugate modes. 
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V. SYSTEMATIC ERRORS 

Systematic errors in measuring the branching fractions mainly originate from the MC 
statistics, the error matrix of the track finding and the 1C kinematic fit, the vertex re- 
quirement, particle identification, the photon efficiency, the a selection criterion, the total 
number of i]){2S) events, and the fitting of the signal. 

1. The MDC tracking efficiency was measured using channels like J/ip — > AA and 
ip{2S) — > 7T + 7T~J/ifj, J/ip — > It is found that the MC simulation agrees with 
data within (1 — 2)% for each charged track. Therefore, 4% is taken as the systematic 
error for events with two charged tracks. 

2. The difference in V z between the positive and negative charged tracks \V+ — V~\ is 
required to be less than 3 cm, this corresponds to an about three-standard-deviation 
requirement. The effect of it is checked with J/i/j — > pnn~ + c.c. candidate events, it is 
found that MC simulates data within 0.6%. This is taken as systematic error of this 
selection criterion. 

3. The photon detection efficiency was studied with different methods using J/tp — > 
7T + 7r~7r° events [14], and the difference between data and MC simulation is about 2% 
for each photon. We take 2% per photon in the analysis. 

4. Particle identification is used in selecting candidate events, and we take 5% as the 
systematic error [11]. 

5. The systematic error from the 1C kinematic fit should be smaller than for the 4C kine- 
matic fit, since there are fewer constraints. Various studies show that the uncertainty 
of the 4C kinematic fit is around 4% [17], so here we conservatively take 4% as the 
error from the 1C kinematic fit. 

6. The uncertainties of background shapes in ip(2S) — > pnir~ , ip(2S) — > pnir + , ip(2S) — > 
pn-R-n + c.c, ip(2S) -> 7Xco,XcO -> pnir~, and ip{2S) -> 7Xc2,Xc2 -> pnn- are 
estimated to be about 3.4%, 4.0%, 12%, 14%, and 25%, respectively, by changing the 
order of the background polynomial and the fitting range. 
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7. The effect of the requirement on the angle between the nearest neutral cluster and the 
missing momentum direction of all charged tracks is checked with J/tp — > pnn~ + c.c. 
candidate events, where the statistics are much higher and the background is much 
lower. We require the same selection criteria on these two channels as ip{2S) — > 
pnix~ + c.c. except for the requirement on a. By applying the requirement on a, we 
measure the efficiency of this selection criterion experimentally. 

The efficiency difference for J/tp — > pnn~ between data and MC simulation is measured 
to be = (101.60 ± 0.53)% . We take f\ = 1.02 as the efficiency correction factor 
for channels containing a n. 

The efficiency difference for J/ip — > pmr + is measured to be = (89.80 ± 0.49)%. 
We take f 2 = 0.90 as the efficiency correction factor for channels containing a n. The 
big difference in the efficiency is due to the fact that the simulation of the hadronic in- 
teraction of the neutron with the detector material is not very reliable. The systematic 
error on the efficiency associated with this requirement is taken as 0.5% . 

8. The energy observed in the BSC associated with the positive charged track is required 
to be less than 0.7 GeV. The effect of this selection criterion is checked with J/ip — > 
pmr + candidate events, the efficiency difference between data and MC simulation is 
measured to be = (100.21 ± 0.27)%, and we take 0.5% as the systematic error on 
the efficiency associated with this selection criterion. 

9. In order to select ip(2S) — > JXcj,Xcj —> pfnT~TT°, ^{2S) — > pfni~7i , and ip(2S) — > 
IXcj, Xcj pfm~ candidate events, four, three, and two photons are required, respec- 
tively. The effect of these selection criteria are checked with J/ip — > pfm~ candidate 
events, and the efficiency difference between data and MC simulation is measured to 
be f^- = (93.8 ± 2.2)%. Because the MC simulates less fake photons than in data. 
We take fs = 0.94 as the efficiency correction factor for three background channels 
analyzed, and take 2.4% as the systematic error on the efficiencies associated with 
these three selection criteria. 

Table I lists all the systematic errors, and the total systematic errors for ip(2S) — > pnn~ , 
ip(2S) -> pnir + , ip(2S) -> pnTT^Tr , i/j(2S) -> 7X c0 , Xco -»• pun', and i^(2S) -> 7X c2 , 
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Xc2 — > P^tt are 8.4%, 8.6%, 16%, 17%, and 27%, respectively. The systematic error for 
ip(2S) -> 7X c j, XcJ -»■ pnn-Ti is 11%. 

TABLE I: Summary of systematic errors (%), the errors common to all modes are only listed once, 
and "• • • " means no contribution. 
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VI. RESULTS AND DISCUSSION 

For the signal channel tp(2S) — > pnn~ + c.c, the branching fractions of the main back- 
ground channels are given above. From the inclusive and exclusive MC simulations, the 
contributions of other background channels in the missing mass distributions of pir~ and 
p7i + show no peak. By fitting the missing mass distributions of pir~ and pir + with signal 
histograms obtained from MC simulation, normalized histograms of the MC simulated main 
background channels measured in this analysis, and first order Legendre functions to de- 
scribe the other backgrounds, 921 ±40 and 914 ±42 events are obtained. The fits are shown 
in Fig. 1. The efficiency for ip{2S) — > pnix" is e 5 = (26.48 ± 0.14)%, and the branching 
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fraction is: 



B(tp(2S) -> pmr - ) 



N sig 
pnn~ 



(2.45 ±0.11 ±0.21) x 10" 



The efficiency for ip(2S) — > pn7r + is e 6 = (28.83 ± 0.15)%, and the branching fraction is: 



B{ip{2S) -> pnvr" 



(2.52 ±0.12 ±0.22) x 10" 



N^{2S) • €6 ■ h 

Here /i,/2 are the correction factors to the efficiencies, and the first errors are statistical 
and the second ones are systematic. 

Taking events with the missing mass within ±0.1 GeV/c 2 around the neutron mass, we 
get 851 and 849 for ip(2S) — ■> pnn^ and i/)(2S) — > pn7r + events, respectively. The Dalitz plots 
for these two channels, which are similar, are shown in Fig. 6. The asymmetry between pn 
and nil is partly due to the difference in detection efficiency and may be partly due to isospin 
symmetry breaking effects from the electromagnetic interaction [18]. From the plots, the 
contribution of the N* states at around 1.4 — 1.5 GeV/c 2 can be seen, and there is a possible 
vertical band at around m 2 = 4.75 (GeV/c 2 ) 2 . A similar band in the horizontal direction is 
even less clear since the events on the right side, where the recoil proton (antiproton) has 
low momentum and cannot be detected well by the detector, have low efficiency. 
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FIG. 6: Dalitz plots for (a) ^(25) -► putt' and (b) i/)(2S) 
0.938| < 0.1 GeV/c 2 . 



prnr 
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Figure 7 shows the pn~ (or pn + ) and nn~ (or nir + ) invariant mass distributions for 
ip(2S) — > pn7Y~ + c.c. In order to investigate the behavior of the amplitude squared as a 
function of invariant mass, the invariant mass distributions are divided by phase space and 
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corrected by the mass dependent efficiency. The results are shown in Fig. 8. There is a 
large accumulation of events below 1.5 GeV/c 2 , which may be due to iV*(1440), iV*(1520), 
iV*(1535), etc. The cluster of events above 2 GeV/c 2 is partly due to the reflection of the 
iV*(1440) etc., and partly may due to high mass N* states, for example TV* (2190), iV*(2220), 
TV* (2250). No clear iV*(2065) peak is observed in the plots, although we can not rule out 
its existence. Partial wave analysis is necessary to obtain more information about the N* 
components in the data. However, here the statistics are not large enough for such an 
analysis. 
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FIG. 7: pir (or pir + ) and mr (or nir + ) invariant mass distributions for ip(2S) — > pnir + c.c. 
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FIG. 8: Data corrected by MC simulated efficiency and phase space versus pn (or pir + ) and nir 
(or nir + ) invariant mass for tp(2S) — ► pn7r~ + c.c. candidate events. 
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Using the J ftp decay results in Refs. [19-21], we obtain: 



^- = B B<fn ^"? = (12.0 ±i.5)%, 

B(J/ip — > pun-) 

which agree with the 12% rule within one standard deviation. 

Using the branching fraction of ip{2S) — > pp7r° from Ref. [12], the ratio of B{jp{2S) — > 
ppir ) : B(if>(2S) -> prar~) : B(if>(2S) -> pn7r+) is measured to be 1 : (1.86 ± 0.27) : 
(1.91 ± 0.27). In calculating this ratio, the common systematic errors between them have 
been removed. This ratio is consistent with the ratio 1:2:2 predicted by isospin symmetry. 



VII. SUMMARY 

Using 14 million ip(2S) events, the branching fractions of ip{2S) — > pnn~ , ip(2S) — > 
pn7T+, ^(25) -> pnn-ir°, ip(2S) -> 7X c0 , XcO -> pmr - , and ^(25") -> 7Xc2,Xc2 -> P^vr" are 
measured to be: 

B{i){2S) -> pnTT") = (2.45 ± 0.11 ± 0.21) x 10~ 4 , 

B(ip(2S) -> pnvr + ) = (2.52 ± 0.12 ± 0.22) x 10" 4 , 

B{i){2S) -> pnvr-Tr ) = (3.18 ± 0.50 ± 0.50) x 10" 4 , 

5(^(25) -> 7Xco, Xco - pnn~) = (1.10 ± 0.24 ± 0.18) x 10" 4 , 

B(ijj(2S) -> 7X C 2, Xc2 -> P^O = (0.97 ± 0.20 ± 0.26) x 10" 4 , 

where the first errors are statistical and the second are systematic. The upper limit of 
J2 2 J=0 B(i/j(2S) -> 7x c j, XcJ -> puTT'ir ) is estimated to be 1.2 x 10" 4 at the 90% C. L. 
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